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Abstract 


Based on efe annihilation data samples collected with the BESIII detector at the BEPCII collider at 13 center of 
mass energies from 4.008 to 4.600 GeV, measurements of the Born cross section of e* e^ — ppn? are performed. 
No significant resonant structure is observed in the measured energy-dependent cross section. The upper limit on the 
Born cross section of e*e^ — Y (4260) — ppn? at the 90% C.L. is determined to be 0.01 pb. 

Keywords: hadrons, cross section measurements, Y (4260) 


1. Introduction 


The Born cross section of ete~ — ppr” in the 
vicinity of the (3770) has been measured recently 
by BESIII [1]. Information on the cross section of 
e*e- — ppn? at higher energies is however still lack- 
ing. The experimental data on the cross section of 
e*e- — hadrons can be used as an input to calcu- 
late the hadronic vacuum polarization via dispersion in- 
tegrals [2, 3, 4, 5]. 

The charmonium-like state Y (4260) was first ob- 
served in its decay to «^ J/w [6]. So far, there 
is no evidence of the Y (4260) in the measured open 
charm decay channels [7, 8] and R value scans [9, 10, 
11, 12, 13, 14, 15]. Many theoretical models have 
been proposed to interpret the nature of Y (4260), e.g. 
as a tetraquark state [16], a DD or DoD* hadronic 
molecule [17], a hybrid charmonium [18, 19], or a bary- 
onium state [20]. Searches for new decay modes of 
the Y (4260) may provide information that can shed 
light on the nature of Y (4260). In particular, the hy- 
brid model [18] predicts a sizable coupling between the 
Y (4260) and charmless decays. 

In this analysis, we report measurements of the cross 
section of efe — ppr? based on the e*^e- annihi- 
lation samples collected with the BESIII detector at 13 
center of mass energies in the range /s = 4.008— 4.600 
GeV as shown in Table 1. Results of the measure- 
ments can be used to estimate the cross section of pp — 
Xeen, which is of high importance for the planned 
PANDA experiment [21] at FAIR in Darmstadt, Ger- 
many. 
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2. BESIII detector and Monte-Carlo simulation 


The BESIII detector [22] is a magnetic spectrometer 
operating at BEPCII, a double-ring efe collider with 
center-of-mass energies between 2.0 and 4.6 GeV and 
a peak luminosity of 1033 cm ?s^! near the (3770) 
mass. The cylindrical core of the BESIII detector con- 
sists of a helium-based main drift chamber (MDC), a 
plastic scintillator time-of-flight system (TOF), and a 
CsI(TI) electromagnetic calorimeter (EMC) that are all 
enclosed in a superconducting solenoidal magnet pro- 
viding a 1.0 T magnetic field. The solenoid is sup- 
ported by an octagonal flux-return yoke with resistive 
plate counter muon identifier modules interleaved with 
steel. The acceptance for charged particles and photons 
is 93% of the 47 solid angle, and the charged-particle 
momentum resolution is 0.596 for transverse momenta 
of 1 GeV/c. The energy resolution for showers in the 
EMC is 2.5 (596) for 1 GeV photons in the barrel (end- 
caps) region. 

A GEANT4-based [23] Monte Carlo (MC) simulation 
software package is used to optimize the event selec- 
tion criteria, estimate backgrounds and determine the 
detection efficiency. For each energy point, we generate 
200,000 signal MC events of ete~ — ppr? uniformly 
in phase space. Effects of initial state radiation (ISR) 
are simulated with KKMC [24], where the line shape of 
the production cross section of e* e^ — ppr is taken 
from results of the measured cross section iteratively. 
Effects of final state radiation off charged particles are 
simulated with PHOTOS [25]. 
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To study possible backgrounds, a MC sample of in- 
clusive Y (4260) decays, equivalent to an integrated lu- 
minosity of 825.6 pb~!, is also generated at \/s = 4.26 
GeV. In these simulations, the Y (4260) is allowed to 
decay generically, with the main known decay chan- 
nels being generated using EVTGEN [26] with branch- 
ing fractions set to world average values [27]. The re- 
maining events associated with charmonium decays are 
generated with LUNDCHARM [28], while the continuum 
hadronic events are generated with PYTHIA [29]. QED 
events (e*e^ — ete, tp, and yy) are generated 
with KKMC [24]. The sources of backgrounds at other 
energy points are assumed to be similar. 


3. Event selection 


The final state in this decay is characterized by two 
charged tracks and two photons. Two charged tracks 
with opposite charge are required. Each track is re- 
quired to have its point of closest approach to the beam 
axis within 10 cm of the interaction point in the beam 
direction and within | cm of the beam axis in the plane 
perpendicular to the beam. The polar angle of the track 
is required to be within the region of | cos 0| < 0.93. 

The time-of-flight and the specific energy loss 
dE / dx of a particle measured in the MDC are combined 
to calculate particle identification probabilities for pion, 
kaon, and proton hypotheses. For each track, the parti- 
cle type yielding the largest probability is assigned. In 
this analysis, one charged track is required to be identi- 
fied as a proton and the other one as an anti-proton. 

Photon candidates are reconstructed using clusters 
of energy deposited in the EMC. The energy deposited 
in nearby TOF counters is included in EMC measure- 
ments to improve the reconstruction efficiency and the 
energy resolution. Photon candidates are selected by re- 
quiring a minimum energy deposition of 25 MeV in the 
barrel EMC (| cos 0| < 0.8) or 50 MeV in the end cap 
EMC (0.86 < |cos0| < 0.92). To reject photons radi- 
ated from charged particles, the angle between the pho- 
ton candidate and the proton is required to be greater 
than 10 degrees. A more stringent cut of 30 degrees 
between the photon candidate and the anti-proton is ap- 
plied to exclude the large number of photons from anti- 
proton annihilation. 

For events with one proton, one anti-proton, and at 
least two photons, a kinematic fit (4C) with the total four 
momenta of all particles constrained to the energy and 
three momentum-components of the initial e* e^ sys- 
tem is applied. When more than two photons are found 
in an event, all possible ppyy combinations are consid- 
ered and the one yielding the smallest x26 is retained 


for further analysis. The Xic is required to be less than 
30. After selecting the ppyy candidate, the 7° candi- 
dates are selected by requiring |M (yy) — m,o| < 15 
MeV/c?, where m,» is the nominal 7° mass [27]. 

The Dalitz plot for the events passing the above se- 
lection criteria for data at \/s = 4.258 GeV is shown 
in Fig. 1(a). The corresponding invariant mass spectra 
of pp, px? and pz? are shown in Fig. 1(b), (c) and (d), 
respectively. 

The potential backgrounds for e*e^ — ppr” are 
studied using the inclusive MC sample at ys = 4.26 
GeV. After imposing all event selection requirements, 
the remaining background events are found to have the 
final state topologies efe — pp, yypp and yyypp. 
No other background survives. The non-7° background 
events can be evaluated from events in the 7° side- 
bands. The 7° sideband regions are defined as 0.07 < 
M(yy) < 0.10 GeV/c? and 0.17 < M(yy) < 0.20 
GeV/c?. The background contamination estimated us- 
ing 7° sidebands at v/s = 4.258 GeV is 0.3%. The back- 
ground contributions are neglected in the subsequent 
analysis. 
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Figure 1: (a) Dalitz plot for the selected et e^ — ppr? candidates of 
data and invariant mass spectra of (b) pp, (c) px? and (d) pz? at Js 
= 4.258 GeV. In (b), (c) and (d), the points with error bars show data 
and the red histograms show MC projections of partial wave analysis 
fit described in the text. 


4. Study of intermediate structures by Partial Wave 
Analysis 


As shown in Fig. 1, a prominent structure near the 
threshold in the pp mass spectrum is visible. Structures 
are also seen in the pz? and pz? mass spectra. To evalu- 
ate the detection efficiencies of the decay e^ e^ — ppr? 
properly, a partial wave analysis (PWA) is performed 


with the e*e- — ppr” candidates to study the inter- 
mediate states present. 

For the process efe — ppr®, the isospin of the 
ppr? system can be J = 0 or I = 1. The quasi- 
two-body decay amplitudes in the sequential decay pro- 
cesses ete” — pN*(pN*), N*(N*) — pn" (pn^), 
ete — pA*(pA*), A*(A*) — pn? (pr^), e*e- 5 
p* (w*)n9, p*(w*) — pp are constructed in the covari- 
ant tensor formalism [30, 31]. All 177 and 377 states 
above pp threshold, N* and A* states with spin up to 
5/2, listed in the summary tables of the PDG [27], are 
considered in this analysis. According to the framework 
of soft meson theory [32], the off-shell decay process 
should be included. Thus, N (940) with a mass of 940 
MeV/c? and zero width representing a virtual proton 
which could emit a 7° is considered as a possible com- 
ponent. No isoscalar vector meson is considered, since 
there is no candidate above the pp threshold in the sum- 
mary tables of the PDG. The p* states are parameter- 
ized by a constant-width relativistic Breit-Wigner (BW) 
propagator with barrier factors included. The N* and 
A* states are parameterized by a BW propagator as de- 
scribed in Ref. [30]. The resonance parameters are fixed 
according to previous measurements [27] due to limited 
statistics. The complex coefficients of the amplitudes 
are determined by an unbinned maximum likelihood fit. 
The details of the likelihood function construction can 
be found in Ref. [33]. 

For p* states with J = 1, the pp final state inter- 
action (FSI) effect using the Jülich model [34] is taken 
into consideration by factorizing the partial wave am- 
plitude into the amplitude without the FSI effect and the 
S wave pp scattering amplitude in the scattering length 
approximation given in Ref. [34]. The direct process of 
ete — ppr” can be modeled by 177 or 377 phase 
space of the pp system (177 or 3 PHSP). All com- 
binations of the components in Ref. [35] are evaluated. 
The changes in the negative log-likelihood (NLL) and 
the number of free parameters in the fit with and with- 
out a resonance are used to evaluate its statistical sig- 
nificance. Resonances with significance greater than 5c 
are retained in the PWA solution. The selection of PWA 
components is performed at the energy points with the 
high statistics, i.e. at \/s = 4.008, 4.226, 4.258 and 
4.416 GeV, as shown in Table 1. The selected com- 
ponents are used to describe the data at other nearby 
energy points. The data at \/s = 4.189 — 4.600 GeV 
can be described by the N (1440), p(2150), p3(1990) 
and 1 PHSP amplitudes. The data at \/s = 4.008 — 
4.085 GeV can be described by the N (1520), N (2570), 
p(2150), p3(1990) and 1 PHSP amplitudes. The 
N (940) is not included in the fits since its significance 


is less than 5c. If we perform an alternative PWA fit 
with N(1440), p(2150), 03(1990) and 1 PHSP at 
s/s = 4.008 GeV, the NLL worsens by 37.8. The change 
of efficiency determined with the alternative fit with re- 
spect to the nominal value is considered as a source of 
systematic uncertainty. Comparisons of the data and the 
fit projection (weighted by MC efficiencies) in terms of 
the invariant mass spectra of pp, pr? and pz? at \/s = 
4.258 GeV are shown in Fig. 1(b), (c) and (d), respec- 
tively. The x? over the number of bins is displayed in 
those figures. 


5. Cross section for e* e^ — ppr? 


The Born cross section for efe — pjpn? is deter- 

mined as 
obs 
B MEME NN (1) 
Co (14-07) (1+8). €: Bro 

where N*'5 is the number of observed events; £ is the 
integrated luminosity; € is the detection efficiency de- 
rived from MC events generated according to the results 
of the PWA fit; (1 + 6”) is the radiative correction fac- 
tor, which is taken from a QED calculation taking the 
line shape of the cross section efe — ppr? of data 
as input in an iterative procedure; (1 + 6”) is the vac- 
uum polarization factor, including leptonic and hadronic 
contributions, taken from a QED calculation with an ac- 
curacy of 0.5% [36]; and Bo is the branching fraction 
of 7° decaying to yy according to the PDG [27]. The 
measured Born cross section of ete~ — ppr’ at each 
energy point is listed in Table 1. 

Uncorrelated systematic uncertainties in the Born 
Cross section measurements mainly originate from the 
1? mass window requirement, kinematic fit and the in- 
termediate states in PWA. The systematic uncertainty 
from the requirement on the 7° signal region is esti- 
mated by smearing the invariant mass of the yy pair 
in the signal MC with a Gaussian function to compen- 
sate for the resolution difference between data and MC. 
The parameters for smearing are determined by fitting 
the 7° distribution of data with the MC shape convo- 
luted with a Gaussian function. The difference in the 
detection efficiency between signal MC samples with 
and without the extra smearing is taken as the systematic 
uncertainty. The systematic uncertainty due to the kine- 
matic fit is estimated by correcting the helix parameters 
of charged tracks for the signal MC sample according 
the method described in Ref. [37]. The difference in the 
detection efficiency between the MC samples with and 
without this correction is taken as the systematic uncer- 
tainty. The systematic uncertainty from the intermediate 


Table 1: The results on e*e- — ppr?. Shown in the table are the integrated luminosity £, the radiative correction factor (1 + ó”), the vacuum 
polarization factor (1 + 6”), the number of observed events N°S, the detection efficiency e and the Born cross section o? (e* e^ — ppr?) at 
each energy point. The errors of e are from the PWA fit. The first errors of o? are statistical, and the second ones are systematic. 


Vs(GeV) £[p | (1-5) (1-6) Ne e[96] a” [pb] 
4.008 482.0 0.967 1.044 1074433 43.9409 5.0943: 0.1870 75 
4.085 52.6 0.992 1.052 106411 43.7414 447-0467?! 
4.189 43.1 1.025 1.056 7529 44A7+10 3644+0.1431019 
4.208 54.6 1.031 1.057 93-10 449+16 3.52+0.1239105 
4.217 54.1 1.034 1.057 82-10 43.4413 3.24: 0.37 + 0.18 
4.226 1047.3 1.037 1.056 1611441 45.2+0.5 3.15 0.08 + 0.14 
4.242 55.6 1.042 1.056 89-9 446-11 3302 0.367019 
4.258 825.6 1.048 1.054 1203435 434+40.5 3.08 + 0.10*0]4 
4.308 44.9 1.063 1.053 538258 — 46.0:E14 . 2.924:0.3370 15 
4.358 539.8 1.081 1.051 668426 447-11 24840.1111 
4.387 55.2 1.087 1.051 5728 475-18 1.92 0.26 0.10 
4.416 1028.9 1.098 1.053 1133434 446406 2.62 0.10*019 
4.600 566.9 1.124 1.055 — 474422 43.8: 0.8 1.63 4 0.08 + 0.08 


states in PWA includes those from the BW parametriza- 
tion, resonance parameters and extra resonances. Un- 
certainties from the BW parametrization of intermedi- 
ate states are estimated by replacing the BW formula 
of N(1440) and N(1520) as used in Ref. [30] with a 
constant BW formula and replacing those of (2150) 
and p3(1990) with the BW formula with the Gounaris- 
Sakurai (GS) model [38]. In the PWA fit, the resonance 
parameters are fixed according to the previous measure- 
ments [39, 40]. Alternative fits are performed in which 
the resonance parameters are set as free parameters and 
the changes in the results are taken as systematic uncer- 
tainties. Uncertainties from additional resonances are 
estimated by adding the most significant additional res- 
onance among each J? assignment in Ref. [35] into the 
PWA solution individually, and their influences on the 
cross section measurements are taken as the systematic 
uncertainties. 


Correlated systematic uncertainties among the dif- 
ferent energy points include those from luminosity mea- 
surement (1.0%) [41], MDC tracking (2% for two 
charged tracks) [42], particle identification (2% in total 
for proton and anti-proton) [43], photon detection ef- 
ficiency (296) [44] and radiative correction. The differ- 
ence in e(14-ó") between the third and fourth iteration is 
taken as the systematic uncertainty due to the radiative 
correction, as the radiative correction dependent quan- 
tity «(1 + ó”) converges after three iterations. 


The total systematic uncertainty of the different en- 
ergy points is calculated by adding the individual uncer- 
tainties in quadrature as shown in Table 2. 


6. Upper limit on c(e*e- — Y (4260) — ppr?) 


Figure 2 shows the measured Born cross section of 
ete — ppr° in the energy region studied in this work. 
No significant resonant structure is observed. The upper 
limit on the Born cross section of efe — Y (4260) > 
ppr’ is determined by a least squares fit of 


2 
xp(i$)| , 
(2) 


to the calculated cross sections. In Eq. (2), ocon and 
cy represent the continuum cross section and resonant 
cross section, respectively, and dcon can be described 
by a function of s, Ocon = C/ s^, where the exponent 
À is a priori unknown. The parameter ¢ describes the 
phase between resonant and continuum production am- 
plitudes. The mass m and width T of the Y (4260) are 
fixed to the PDG values [27]. The values of C, A, oy, 
and the interference phase ¢ are free in the fit. The 
uncorrelated systematic uncertainties in the Born cross 
section measurements are directly considered in the fit 
and the effect of the correlated systematic uncertain- 
ties on the final results is estimated by the method in 
Ref. [45], in which the error propagation is determined 
from shifting the data by the aforementioned correlated 
uncertainties and adding the deviations in quadrature. In 
addition, the uncertainties for the beam energy measure- 
ments of all the data points taken from Ref. [46] are con- 
sidered in the fit. The best fit function is shown in Fig. 2 
as the solid line. The dashed line represents the fit with 
oy = 0. The optimal value of oy is (1.63-5.9) x 10^? pb 
with a statistical significance of 0.50. The significance 
is calculated based on the changes in the x? value and 


r 
a(s) = |V/Ocon + y dn di 


—————— = 
— m? + imr 


Table 2: Summary of systematic uncertainties on the Born cross section of e^ e^ — ppr? (96). 


Sources /\/s (GeV) 4.008 4.085 4.189 4.208 4.217 4.226 4.242 4.258 4.308 4.358 4.387 4.416 4.600 
Luminosity 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 
MDC tracking 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 
PID 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 
Photon detection 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 
Kinematic fit 2.1 1.9 1.8 1.6 1.7 1.6 1.6 2.1 15 1.8 1.5 1.8 1.6 
1? mass resolution 0.2 0.2 0.2 0.2 0.2 0.3 0.2 0.4 0.6 0.4 0.3 0.4 0.4 
Radiative correction 1.9 1.9 1.9 1.9 1.9 1.9 1.9 1.9 1.9 1.9 1.9 1.9 1.9 
Intermediate statesin PWA *13.— *15.— "55 taa — "35 — "13 tio — "$9 — *gg — "32 — "yg — 793 — £19 
Total Dss az —-52  -62  -56 +45 "ig -so -as -4s  -49  -49 +48 
the number of free parameters in the fit with and without = 
the assumption of existence of the Y (4260) resonance. S BL —+— data 
The result for the phase between resonant and contin- S "S continuum 
. ç , [z4 L resonance + continuum 
uum production amplitudes is ç = 3.4+1.0. The param- e 45 
eters describing the slope of the continuum cross section + [ 
are C = (5.4 + 5.3) - 10° GeV?^pb and À = 4.2 + 0.4. fas 3- 
The upper limit on oy at the 90% C.L., oy”, is deter- M r 
E oy? oo a |` 
mined by f” G(oy, ooy )dz/ f Gloy, coy )dz = o 2 
0.9, where G(ocy,osíy) is a Gaussian function with | d pedem te ig 
mean value oy = 1.6 x 107? pb and standard deviation 4 4.2 4.4 4.6 
Goy = 5.9 x 107? pb. The uncertainties from mass and (s (GeV) 


width of the Y (4260) are considered by varying them 
by one standard deviation according to the PDG val- 
ues [27] and the most conservative oy is taken as the 
final result. The obtained upper limit is 0.01 pb. 


7. Summary 


Based on 13 data samples between 4/s = 4.008 and 
4.600 GeV collected with the BESIII detector, the pro- 
cess e*e- — ppr’ is studied. The Born cross section 
of e*e- — ppr? is measured. No resonant structure 
is observed in the shape of the cross section. The upper 
limit on the Born cross section of e*e^ — Y (4260) > 
ppr” at the 90% C.L is estimated to be 0.01 pb. 


Figure 2: Fit to o(e* e- — ppr?) with resonance and continuum 
(solid line), or only continuum term (dashed line). Dots with error 
bars are the measured Born cross sections. The uncertainties are sta- 
tistical only. 
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